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Section  1 
INTRODUCTION 


This  is  the  final  report  by  McDonnell  Douglas  Astronautics  Company  (MDAC)  of 
work  conducted  for  the  Naval  Research  Laboratories;  Washington,  D.C.  under 
contract  N001 73-78-01 61 ,  titled  "Application  of  Nonlinear  Crystals."  An 
oral  presentation  was  made  in  Washington,  D.C.  on  July  9,  1979  and  the 
view  foils  used  are  included  as  Appendix  A-I. 

The  body  of  this  report  is  a  narrative  description  of  the  July  9  presentation 
referencing  Appendix  A-I  figures  and  comments.  The  outline  used  for  the  oral 
presentation  is  given  on  page  A-2.  Auxiliary  theoretical  material  has  been 
included  in  Appendixes  A- I I .  The  topics  are  keyed  into  this  narrative  in 
Section  1.1  and  a  summary  of  the  report  is  given  in  Section  1.2. 

1 . 1  Out! ine 

A  program  overview  including  the  work  statement  and  contract  history  is 
given  in  Section  2.  The  body  of  the  report  starts  in  Section  3  where  prior 
work  on  frequency  conversion  done  primarily  by  the  Russians  is  discussed 
briefly.  In  Section  4  a  detailed  description  of  experimental  work  carried 
out  under  this  contract  is  divided  into  ten  (10)  subsections.  In  Section  4.1 
the  joint  transform  correlator  work  with  HeNe  is  described.  This  is  followed 
in  4.2  by  a  description  of  MDAC's  newly  acquired  frequency  conversion  facilities. 
Section  4.3  is  an  analysis  of  the  optical  properties  of  NRL  bit  shifted  photo¬ 
graphic  transparencies  supplied  for  use  on  this  project.  A  variety  of  experi¬ 
mental  factors  which  contribute  to  upconversion  efficiency,  and  resolution  are 
described  in  Section  4.4  to  4.9  in  the  framework  of  upconversion/corvol ution 
operations.  Last,  in  Section  4.10,  an  experimental  test  called  the  N-pattern 
test  is  described  and  i llustrated  for  the  evaluation  of  second  harmonic  generation 
(SHG)  conditions  as  a  function  of  temperature,  angle,  and  SHG  crystal  translation 
and  rotation.  Next,  upconversion  and  downconversion  conditions  for  real  time 
convolution  and  correlation  operations  are  discussed  in  Section  5.  Finally, 
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the  conclusions  on  applications  of  nonlinear  crystals  for  up-  and  downconversion 
are  given  in  Section  6.1  and  recommendations  for  further  study  are  given  in 
Section  6.2.  Appendix  A-I  is  the  collection  of  view  foils  used  for  the  NRL 
Program  Review  in  Washington,  D.C.  Appendix  A-II  is  a  review  of  the  basic 
theory  behind  frequency  conversion,  the  correlation  operator,  and  the 
correlation  of  optical  data  using  nonlinear  crystals. 

1 . 2  Summary 

Although  personnel  changes  and  experimental  problems  delayed  the  program, 
the  main  goals  which  were  reached  included  establishing  a  diffraction  limited, 
Nd:YAG  facility  for  investigating  nonlinear  crystal  interactions  for  applica¬ 
tions  to  ODP,  developed  experimental  procedures  for  optimizing  ODP  optical 
conversion  efficiency,  and  defined  the  experimental  conditions  for  performing 
two  frequency  correlation  by  both  upconversion  and  downconversion.  Equip¬ 
ment  for  three  wavelength  upconversion  correlation/convolution  are  available 
for  preliminary  tests  of  this  more  complicated  ODP  process.  Finally,  we  are 
ready  to  process  the  NRL  bit  shifted  data  to  measure  correlation  accuracy  and 
Space-Bandwidth  Product  (SBWP). 

An  N-pattern  test  was  developed  to  determine  optimum  temperatures  for  using 
CD*A  as  a  Fourier  plane  nonlinear  crystal  in  up  and  downconversion.  By 
measuring  correlation  peak  widths  using  N-pattern  test  data  the  system 
resolution  was  found  to  be  11  tp/mm  which  was  the  diffraction  limit  of  the 
lens.  A  50  mm,  f/4  lens  is  available  to  substitute  in  the  450  mm  lens  now 
in  use.  This  will  increase  the  resolution  by  a  factor  of  9  to  100  tp/mm. 

Correlation  and  convolution  can  be  interchanged  by  rotating  real  patterns. 

This  is  done  by  rotating  180°  one  of  a  pair  of  input  patterns  which  is  equiva¬ 
lent  to  complex  conjugation  in  the  Fourier  plane  domain.  A  test  complex 
conjugate  input  pattern  was  prepared  and  is  now  ready  to  use  for  correlation 
by  upconversion. 

By  using  filters  at  the  input  and  image  planes,  it  is  possible  to  identify 
the  ODP  terms  and  the  different  components  contributing  to  signal  and  noise. 

In  this  way  we  were  able  to  find  the  conditions  for  phase  matching  by  tempera¬ 
ture  and  angle  tuning  for  downconversion  using  the  N-pattern  test. 


Section  2 
PROGRAM  OVERVIEW 


2. 1  Statement  of  Work 

A  task  breakout  of  topics  to  be  covered  is  given  on  page  A-5.  In  Task  I 
correlation  of  two  images  was  proposed  using  convolutional  upconversion  of 
Hermitian  functions  which  in  this  special  case  is  equivalent  to  the  correla¬ 
tion.  The  Space-bandwidth  product  (SBWP)  would  then  be  measured  from  the 
correlation  peak  width  at  the  output.  This  SBWP  would  then  be  compared, 
in  Task  II,  with  comparable  SBWP  values  determined  using  an  optical  correla¬ 
tor  and  digital  computer  techniques.  In  Task  III,  the  correlation  which  is 
generated  as  a  bulk  zero  frequency  field  would  be  read  out  using  a  third 
frequency  (color)  HeNe.  Task  4  was  designed  to  demonstrate  correlation  of 
NRL  codes  in  real  time,  find  the  bit  shift  accuracy  by  direct  measurement 
and  determine  maximum  bit  density  limitations.  Last,  Task  5  was  to 
select  crystals  with  large  nonlinear  coefficients,  determine  optimum 
pulse  matching  conditions,  develop  equations  to  characterize  field  angle, 
and  then  select,  specify  and  purchase  upconversion  crystals  for  the 
experimental  portion  of  this  program. 


Historical  Summary 


The  program  history  is  summarized  on  page  A-6  and  in  Figure  1.  Briefly, 

Peter  Guilfoyle  left  the  program  in  January,  1979  after  the  upconversion 
and  laser  crystals  were  ordered,  and  he  and  Carol  Kastner  completed  the 
Joint  Transform  Correlator  (JTC)  laboratory  work  in  January.  C.  Kastner 
then  assembled  the  equipment  and  started  upconversion  experiments  alone  until 
May  1979  when  J.  Grossman  joined  the  program.  After  correcting  operational 
laser  problems,  some  equipment  modifications  were  implemented  and  the 
convolution  work  completed. 
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Section  3 

PRIOR  WORK  ON  CORRELATION  AND  CONVOLUTION 


Eremeeva,  et  al.  (Reference  1)  published  results  in  1976  showing  that  optical 
convolution  and  correlation  of  two  scenes,  illuminated  by  two  optical 
frequencies  and  respectively,  could  be  accomplished  using  non-linear 
crystals.  The  convolution  operation  is  obtained  when  the  two  frequencies 
are  the  sum  (  ^  +  u^)  frequency  and  the  correlation  operation  is  obtained 
when  the  two  frequencies  are  the  difference  (u^  -  u~)  frequency.  They 
published  results  (reproduced  in  Figure  A-2),  wherein  they  show  the  results 
of  interaction  of  circular  apertures  illuminated  by  combinations  of  the 
fundamental  color  ^ ,  the  second  harmonic  ^  and  the  third  harmonic  Xy 
Inset  a^  shows  the  three  input  plane  hole  symmetry  patterns  used.  A  com¬ 
panion  background  discussion  is  given  in  Sections  1  to  3,  Appendix  A-2. 

Eremeeva  calls  case  d  auto-convolution  because  the  same  wavelength  X-j  is 
used  to  illuminate  both  holes  which  we  distinguish  as  F  and  G,  the  Fourier 
transforms  of  scenes  f  and  g  at  the  input  plane.  Since  the  images  are  real, 
the  Fourier  tranforms  also  contain  the  complex  conjugate  functions  F*  and 
G*.  Experimentally  this  is  accomplished  as  shown  schematically  in  Figure  3d_. 
The  laser  1_  emits  the  fundamental  X-|  through  a  beam  expander  (in  1_)  and 
illuminates  the  object  plane  4  through  a  filter  3  which  passes  only  X-j .  The 
object  plane  4  and  a  second  harmonic  generating  crystal  6^  are  at  conjugate 
focii  of  the  Fourier  transforming  lens  _5.  Thus,  the  Fourier  transforms 
(F  +  F*)  and  (G  +  G*)  of  the  upper  and  lower  object  plane  light  distributions 
are  multiplied  by  the  non-linear  crystal  interaction  to  produce  terms  FF*, 
GG*,  and  FG*  and  F*G  at  w  ■  0  (the  DC  terms)  and  FF,  GG,  F*F*,  G*G*,  FG 
and  F*G*  at  w=  2co^  =  c/X£  where  c  is  the  velocity  of  light  (Reference  2). 

A  second  Fourier  Transforming  lens  5^  images  the  convolutions  and/or  correla¬ 
tions  (as  the  case  may  be)  on  image  plane  ]_  through  a  second  filter  3_  which 
in,  Figure  3^,  passes  only  X7. 
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It  is  easy  to  identify  the  locations  of  the  different  terms.  The  interaction 
of  the  light  through  the  upper  aperture  F  with  itself  to  produce  the  real 
convolution  FF  +  F*F*  at  X2  appears  at  the  lower  position  in  the  image 
plane.  Similarly,  the  terms  GG  +  G*G*  appear  at  the  upper  image  position. 

The  cross-convolution  FG  +  F*G*  appears  at  a  position  halfway  between  the 
two. 

The  energy  which  is  in  the  correlation  terms  goes  into  a  local  static 
deformation  of  the  crystal  at  the  Fourier  plane  which  then  relaxes  at  a 
rate  controlled  by  potential  energy  dissipation  by  phonon  emission.  It  was 
postulated  that  this  field  might  be  interrogated  by  a  different  frequency 
A..-  shown  in  Figure  3a.  Taking  into  account  the  crystal  distortion  life¬ 
time,  the  stress  should  be  interrogated  with  a  pulse  at  wavelength 
synchronized  with  the  fourier-plane  stress-generating  pulse  at  A-,.  Since 
the  ions  in  the  crystal  cannot  respond  to  the  electromagnetic  generating 
field,  a  time  delay  will  be  required  to  permit  the  electronically  stored 
DC  field  energy  to  transfer  into  lattice  deformation  before  the  DC  field 
can  be  read  out. 

In  correlation  and  convolutions  of  Figures  2b  and  c  and  3£  and  c,  the 
third  harmonic  A3  is  also  used.  In  this  case  a  single  upconversion  crystal 
2  is  used  to  generate  A-|  and  Ao  at  the  input  plane  for  convolution  (case  £) 
and  crystals  2  and  3  are  used  to  generate  A-j  and  A3  at  the  input  plane  for 
correlation  (case  bj. 

The  different  combinations  of  input  terms  F.  and  G-  which  are  connected 

*  J 

with  the  A.j  and  Aj  for  correlation  and  convolution  are  summarized  on  the 
right  hand  side  of  Figure  2.  The  companion  designations  correlation, 
auto-correlation,  convolution  and  auto-convolution  are  the  Russian  workers 
designation  for  the  process  of  using  different  or  the  same  wavelengths  for 
downconversion  (correlation)  or  upconversion  (convolution). 

A  terminology  problem  becomes  evident  when  we  try  to  consider  both  wave¬ 
length  and  scene  convolution/correlation.  The  Russians  imply  wavelength 
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and  common  usage  implies  pattern.  This  is  shown  on  page  A- 1 1 .  There  are 
essentially  four  cases  shown  in  the  two  by  two  matrix.  In  what  follows  we 
propose  to  use  the  following  rules.  Correlation  (convolution)  and  cross¬ 
correlation  (cross-convolution)  are  synonomous.  Auto-(  )  refers  to 
pattern  convolution.  Auto-A-  (  )  refers  to  down  (up-)  conversion 

using  the  same  wavelength.  Finally,  auto-auto-A-  (  )  will  be  called 

A  -(  )  when  a  pattern  is  correlated  (convolved)  with  itself  by  SHG. 


Section  4 
EXPERIMENTAL 


4. 1  Joint  Transform  Correlator  (JTC) 

As  part  of  Task  II  (page  A-5),  23  x  23  bit  arrays  both  zero-shifted  and 
128  bit-shifted  input  patterns  were  correlated  using  the  MDAC-JTC,  Figure  4. 
These  patterns  and  their  output  have  been  published  (Reference  3)  and 
discussed  in  detail  there  (See  also  Section  4,  Appendix  All  for  some  recent 
analytical  work).  These  results  will  be  compared  at  some  future  time  with 
the  equivalent  real  time  data  derived  from  optical  processing  with  CD*A 
crystals. 

4.2  Frequency  Conversion  Facilities 

The  Holobeam  5050Q  Nd:YAG  Laser  was  purchased  based  on  lower  cost,  and 
quickest  delivery.  The  Q-switched  oscillator  amplifier  combination 
generates  65  mj  per  pulse  with  a  20  Hz  repetition  rate. 

CD*A  was  found  to  be  the  best  choice  of  crystals  presently  available  for 
SHG.  Three  different  length  crystals,  5  mm,  10  mm  and  20  mm,  were 
purchased.  Each  having  the  same  cross  section  (15  mm  x  15  mm)  cut  perpen¬ 
dicularly  to  the  optic  axis  for  temperature  tuning.  A  smaller  crystal  of 
RDP  used  for  third  harmonic  generation  (THG)  is  available  for  three 
frequency.  It  is  too  small  and  have  convolution  without  excessive 
vignetting  if  used  for  down  correlation  in  the  Fourier  plane. 

The  experimental  equipment  and  setup  are  shown  in  the  annotated  Figures  5, 

6,  and  7. 

4. 3  NRL  Transparencies 
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4.3.1  Resolution  Tests  with  HeNe  Laser 

NRL  sent  transparencies  of  pseudo-random  codes  to  MDAC  to  be  used  for 
measuring  cross-correlation  parameters  using  non-linear  crystal  interac¬ 
tions.  As  a  first  step,  optical  transmission  and  resolution  tests  were 
made  using  a  HeNe  laser  with  the  CD*A  crystal  in  and  out  of  the  Fourier 
plane.  Some  of  these  test  results  are  shown  in  Figure  8.  Two  points 
are  worth  particular  attention. 

First,  we  find  that  the  binary  black  and  white  imagery  has  a  noise 
component  where  the  black  is  not  quite  black  Figure  8a_,  10  sec  exposure, 
and  the  white  is  not  quite  white,  Figure  8d_  and  e,  1/100  sec  exposure. 

The  existence  of  this  "noise"  component  should  be  included  in  future  S/N 
measurements  of  ideal  resolution. 

Second,  multiple  images  appear,  compare  Figures  8a^  and  8b_,  when  the  CD*A 
crystal  is  placed  in  the  Fourier  plane.  The  most  likely  source  is  multiple 
reflection  between  the  crystal  surfaces  and  the  input  plane  transparency. 

This  one  and  others  are  important  sources  of  noise  which  must  be  minimized 
or  eliminated  for  optimum  system  operation  (see  Figures  11  and  12  following). 

4.3.2  Convolution  SHG  and  Transparency  Rotation  Tests 

The  A  convolutions  (Types  FF  +  F*F*  and  GG  +  G*G*)  are  found  to  be  real 
because  they  are  independent  of  transparency  rotation  (Figures  9a,  b, 
and  c).  However,  the  wander  of  the  faint  pattern  between  the  A  -convolutions 
with  respect  to  the  center  of  rotation  shows  it  is  probably  not  the  auto-2- 
convolution  (FG  +  F*G*).  This  is  verified  by  masking  out  alternately  the 
F  and  G  patterns,  Figure  9d  and  e,  noting  the  pattern  remains.  Although  a 
part  of  the  FG  pattern  wanders,  the  appearance  of  the  edge  cutoff  in  d. 
and  e^  suggests  that  the  most  likely  source  is  the  partial  grey  level 
transmission  observed  in  Figure  8a,  b,  and  c.  Verification  of  these 
alternatives  were  not  pursued. 
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Two  points  are  worth  emphasizing.  First,  the  effects  of  non-uniform 

2 

object  plane  illumination  is  apparent  from  the  variation  in  the  A  - 
convolution  intensity  distribution  with  rotation.  Secondly,  a  portion 
of  the  beam  undergoes  A  -convolution  (SHG)  of  the  type  FF,  which  reproduces 
the  image.  This  shall  be  compared  with  auto-A-convol ution  of  the  type 
F'F  of  different  parts  of  the  same  pattern.  This  is  another  level  of 
complexity  to  add  to  the  classification  matrix  on  page  A- 1 1 . 

4.4  Convolutions  with  Square  Apertures 

When  the  CD*A  crystal  is  set  at  102°,  the  optimum  SHG  temperature,  and 
5/16"  input  plane  apertures  are  separated  (center  to  center)  by  1", 

Figure  ICte,  £  and  c,  only  A  -convolutions  are  observed.  However,  when  this 
separation  is  reduced  to  1/2",  the  FG  auto-A-convol utions  appear  Figure  10d_. 

2 

A  number  of  additional  features  are  worth  noting.  The  A  -convolution  in 
Figure  ICte  shows  a  bright  cross.  However,  the  cross  in  d_  is  dark.  Even 
the  multiple  reflections  at  top  have  dark  crosses.  The  main  difference 
between  the  squares  is  the  ones  in  d  are  larger  then  those  in  a  in  the 
ratio  a_:d_  A  90°  rotation  should  normally  produce  no  change  for 

4- fold  symmetry.  Yet,  the  lower  squares  in  £  have  a  2-fold  bright  band 
symmetry  along  the  horizontal  diagonals.  Finally  the  satellite  reflections 
in  c  are  found  to  have  a  pair  of  parallel  dark  bands. 

2 

The  upper  left  A  -convolution  in  Figure  10£  is  sufficiently  dark  than  a 
square  array  of  lines  is  visible  with  density  3  Ip/mm.  This  places  a  lower 
limit  on  the  resolution  (see  Section  4.10.3). 

4.5  Use  of  Image  Plane  Filters  -  On  Two  Crystal  Experiments 

The  configuration  employed  in  this  experimental  arrangement  is  similar  to 
that  shown  schematically  in  Figure  3c  where  both  crystals  2_  and  6  are 
CD*A  set  for  90°  oriented  SHG.  In  other  words,  in  the  relation  (A^  ;  A2) 

( A^) ,  A 2  =A-| . 


When  a  red  filter  is  used  in  the  input  plane  Figure  11c  and  d,  we  see  both 

2  2  ~  ~ 

A  -  and  auto-A-convolutions.  The  A  -convolutions  are  much  lower  in 

intensity  as  shown  in  c  using  low  speed  film.  This  is  to  be  contrasted 

with  the  results  already  observed  in  Figure  9  where  only  A  -convolutions 

were  found. 

Leaving  out  the  red  filter.  Figure  lla^,  the  green  in  the  input  passes 

through  the  system  without  change  imaging  the  input  circular  holes  in  the 

photograph.  By  measuring  the  ratio  of  the  intensities  of  the  green  light 
2 

at  the  A  position  with  and  without  the  red  filter  we  can  determine  the 
2 

A  -convolution  conversion  efficiency.  Similarly,  the  auto-A-convolution 
efficiency  can  be  found.  This  will  be  a  valuable  tool  when  evaluating 
system  operation. 

4.6  Internal  Reflections,  SHG  and  Noise 

Internal  reflections  have  already  been  discussed  (Sections  4.3  and  4.4). 
The  points  of  interest  for  Figure  12  and  13  are  summarized  on  pages  A-26, 
and  A-28.  Reflections  are  an  important  source  of  noise  and  will  have  to 
be  minimized  in  system  design.  Possible  solutions  include: 

a.  Non-parallel  crystal  faces. 

b.  Index  matching  fluid. 

c.  Absorption  filters  in  the  SHG  crystal  Housing. 

d.  Anti-reflection  coatings  which  transmit  both  the  fundamental  and 
second  harmonic. 

4.7  Convolution  Non-linearity 
2 

The  A  -convolution  process  is  basically  a  field  squaring  process  so  that 
areas  in  a  scene  which  have  twice  the  transmission  will  convolve  to  4 
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times  the  intensity.  This  enhancement  is  graphically  demonstrated  in 
the  convolutions  shown  in  Figure  14. 


4.8  Convolution  Peak  Intensity  Ratios 

The  variations  of  the  peak  convolution  intensity  ratios  is  understood  to 

depend  on  the  angular  extent  of  each  input  scene,  the  separation  between 

them,  temperature,  the  direction  of  the  optic  axis  and  the  longitudinal 

position  of  the  SHG  crystal  relative  to  the  Fourier  plane.  Figure  15  is 

an  example  of  CD*A  angle  effects.  The  effect  of  a  10°C  rise  in  temperature 

is  shown  in  Figure  16  where  the  central  spot  is  a_  disappears  in  b  and  the 
2 

A  -convolutions  become  enhanced.  Crystal  translation  and  rotation  are  not 
as  sensitive  for  90"  tuning  as  for  angle  phase  matching.  This  is  discussed 
in  Section  4.10.4  in  relation  to  Figure  22. 

4.9  Crystal  Damage 

Moving  the  image  plane  4.5  cm  (M 0%  of  the  focal  length)  does  not  appreciably 
2 

change  the  A  -  and  auto-convolution  patterns  Figure  17.  However,  it  does 
affect  high  spatial  frequency  scattering  associated  with  crystal  damage  from 
laser  pulse  energies  which  exceeded  the  damage  threshold  of  CD*A  at  the 
Fourier  plane.  For  example,  the  high  frequency  components  associated  with 
the  upper  circular  aperture  in  Figure  17b  disappears  in  £  when  the  image 
plane  is  translated  -4.5  cm  towards  the  image  plane  lens.  Simultaneously, 
scattering  by  damage  appears  at  the  lower  right  hand  aperture.  This 
demonstrates  the  insensitivity  of  convolutions  (and  correlations)  to 
Fourier  plane  crystal. damage  which  would  be  expected  from  visual  inspec¬ 
tion  to  degrade  the  image  more  than  observed.  It  suggests  that  a  crystal 
could  be  subjected  to  considerable  damage  without  seriously  modifying  its 
effectiveness  for  up  and  down-conversion. 

4.10  N-Test  of  CD*A 

Second  harmonic  generation  of  laser  radiation  is  maximized  when  the  vector 
sum  of  the  two  upconverting  photon  wave  vectors  ^  equals  the  second  harmonic 
wave  vector  Since  k(e)  =  (2ir/n(0)  a)  S(0),  then  for  temperature  tuned. 
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colmear  SHG  6  =  90  and  the  refractive  indexes  of  the  fundamental  and 
second  harmonic  are  equal.  This  is  the  usual  A  -correlation  (convolution) 
doubling  condition  which  for  CD*A  occurs  at  102°C.  When  cross-correlations 
are  to  be  performed,  the  input  plane  images  are  usually  separated  by  at 
least  the  width  of  the  scene.  Then,  the  condition  for  phase  matching 
becomes  that  the  projections  of  the  sum  of  the  k-j  vectors  onto  lL  shall  be 
equal.  In  this  case,  the  angle  between  the  central  rays  of  the  two  inter¬ 
secting  cones  is  about  5°  to  6°  and  the  phase  mismatch  is  equivalent  to  a 
refractive  index  mismatch  of  .00274.  This  refractive  index  change  can  be 
obtained  by  lowering  the  temperature.  However,  as  the  temperature  is 

lowered  to  achieve  refractive  index  matching  for  the  auto-A-cross- 
2 

correlation,  the  A  -correlation  index  matching  conditions  are  lost. 

4.10.1  The  N-Pattern  Test 

To  determine  rapidly  the  phase  matching  conditions  as  a  function  of 
temperature,  an  N  was  cut  out  of  black  paper  1-1/2"  high  and  1-3/8"  wide. 
The  width  of  the  lines  forming  the  verticlas  and  diagonal  of  the  N  was 
0.1".  An  input  plane  N  was  illuminated  by  both  1.064  pm  and  532  nm 
radiation  (see  Figure  3c,  at  25°C).  The  results  are  shown  in  Figure  18a, 
b  and  £  for  different  exposures.  The  green  outline  of  the  (backwards)  N 
can  be  seen  as  well  at  wishbone  shaped  images  which  are  the  auto-A- 
convolutions  of  the  legs  of  the  N  with  the  diagonal.  The  joint  of  the 
wishbone  is  the  90°  phase  matched  condition  ualfway  between  leg  and 
diagonal  and  the  arms  of  the  wishbone  are  the  equivalent  angle  matching 
conditions  as  the  separation  between  the  images  change.  When  a  red  filter 
is  used  to  absorb  the  green  falling  on  the  image  plane.  Figures  18£  and 
e,  the  green  light  that  remains  in  the  3  legs  of  the  N  are  A  -convolutions 
of  each  leg  on  itself. 

4.10.2  Temperature  and  SHG 

A  series  of  N-tests  were  conducted  for  different  Fourier  plane  crystal 
temperature,  from  25°C  to  100°C  intervals.  The  results  are  shown  in 
Figures  18,  19,  20  and  21  as  the  temperature  is  increased.  The  joint  of 
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the  wishbone  moves  towards  the  apexes  of  the  N  and  at  100°C,  which  is 
close  to  "ideal"  SHG  temperature,  the  joint  becomes  coincident  with  the 
apex.  The  wishbone  arms  also  come  closer  together  and  move  towards  the 
apex. 

4.10.3  Effect  of  Exposure 

Two  effects  of  exposure  on  the  film  is  seen.  First  the  lines  widen  to 
a  maximum  limit  and  second  additional  segments  of  the  legs  of  the  N  appear 
and  second  order  intensity  wishbones  appear.  In  Figures  18,  19,  and  20 
additional  wishbones  appear  in  between  the  arms  of  the  main  wishbone, 
whereas  in  Figures  20  and  21,  they  appear  on  the  side  away  from  the  apexes 
of  the  N.  This  shows  that  higher  order  terms  in  the  polarization  tensor 
product  with  the  electric  field  are  present  and  might  contribute 
adversely  to  S/N. 

If  the  low  exposure  frames  in  which  the  red  filter  was  used  are  examined 

carefully.  Figures  18d,  19£  and  20d,  a  thin  bright  line  (the  central 
2 

A  -correlation  peak)  approximately  100  urn  wide  is  found  centered  in  the 
legs  and  diagonal  of  the  N.  This  width  is  found  to  be  in  surprisingly  good 
agreement  with  the  expected  resolution  calculated  from  the  system  para¬ 
meters  x/2Af  =  (10  mm/2  x  1.06  x  10~3  x  400  =)  11.8  tp/mm.  By  using  a 
shorter  (100. rrm)  focal  length  Fourier  plane  crystal  lens,  it  should  be 
possible  to  increase  the  resolution  to  40  £p/mm. 

4.10.4  CD*A  Angle  arjd  Translation  Effects 

When  the  CD*A  crystal  was  translated  vertically  3  mm  off  the  optical  axis, 
the  convolution/correlation  images  disappear,  Figure  22c.  This  sharp 
intensity  decrease  starting  at  2.8  mm  can  be  compared  to  crystal  vignetting 
for  the  first  2  mm  of  motion  Figure  22a  and  b_.  The  effect  of  rotating  the 
optical  axis  0.8°  is  seen  to  be  (the  90°  temperature  tuned  wishbone  joint) 
an  intensity  symmetry  change  in  the  wishbone  joint  from  Figure  22b  to  d. 


Section  5 
DISCUSSION 


The  h-pattern  test  showed  that  as  two  input  scenes  are  moved  apart,  the 
optimum  temperature  for  90°  phase  matching  decreases.  Examples  of  this 
are  shown  in  Figure  23  where  at  25°C  using  CD*A  the  A  -convolutions  are 
missing  in  all  cases,  and  the  auto-X-convolution  is  found  in  e  and  f  but 
absent  in  d_.  From  Figure  18,  the  minimum  distance  for  which  cross¬ 
convolution  should  be  observed  is  2.6  mm  +  1  mm.  Since  the  diagonal  of 

2 

the  largest  pattern,  f,  is  only  21  mm,  the  absence  of  the  A  -convolution 
is  not  surprising.  The  cross  convolution  is  also  absent  in  d  since  the 
maximum  width  of  the  NRL  data  is  30  mm.  Consequently,  only  4  mm  from  the 
outside  of  the  two  patterns  are  used  in  the  auto-X-convolution. 

The  interaction  area  can  be  increased  by  increasing  the  temperature  until 

the  minimum  distance  for  color  conversion  is  10  mm.  However,  this  would 

cause  A  -convolutions  since  the  diagonal  is  cm.  Suppression  of  the  A  ■ 

convolution,  keeping  the  center  to  center  pattern  distance  of  20  mm 

constant,  is  possible  if  the  pattern  width  is  decreased  by  [^2  -  l)mm. 

Then  if  the  DC  field  is  read  out  using  another  wavelength  laser,  the 
2 

A  -correlations  will  not  interfere. 

When  one  of  the  input  plane  patterns  is  rotated  180°  (forming  the  complex 
conjugate)  the  convolutions  become  equivalent  to  correlations.  This  was 
found  to  be  true  by  examining  the  A  -convolutions  of  the  N-patterns.  In 
them,  there  is  a  central  white  line  about  100  urn  wide  which  for  the  lens 
size  used  would  give  the  expected  diffraction  pattern  resolution  limit  of 
11  Jtp/mm.  Therefore,  improved  resolution  is  possible  using  a  lens  with 
a  shorter  focal  length.  A  factor  of  9  improvement  in  resolution  18  is 
theoretically  possible  by  using  a  50  mm  lens  at  f  =  1.6. 


Using  uoconversion  there  are  two  ways  to  obtain  correlations.  These  are: 

a.  Temperature  tune  -  one  central  green  pattern  -  FG*  +  F*G 

b.  Angle  tune  -  two  green  patterns  -  FG*  and  F*G 

This  latter  case  is  one  in  which  the  arms  of  the  wishbone  pattern  are 

used  (resolution  of  narrower  angle  tuning  permitting)  one  being  FG*  and 

the  others  F*G.  In  the  former  case  the  1st  order  diffraction  pattern  would 

be  used  to  separate  FG*  and  F*G  from  the  central  image.  In  this  case  the 

2 

suppression  of  the  A  -correlations  discussed  above  is  mandatory. 


Downconversion  interchanges  the  roles  of  correlation  and  convolution. 
Gormal  patterns  (both  oriented  the  same  direction)  produce  correlation  and 
rotated  oatterns  produce  convolution.  When  only  two  wavelengths  are  used 
(the  fundamental,  Xp  and  second  harmonic, X2)  with  normal  patterns, 
correlations  are  produced  (x-|,  X-j  )-*-  (DC) 

c.  Temperature  tune  -  DC  pattern 

d.  Angle  tune  -  DC  pattern 

The  readout  is  affected  by  a  second  laser,  at  any  convenient  wavelencth, 
wnich  is  synchronized  with  and  delayed  in  time  from  the  IP,  laser  drive 
pulse  to  allow  maximal  transfer  of  stored  slectron  energy  into  a  stress 
field  which  scatters  the  readout  light  pulse. 


Another  correlation  process  is  assymetric  downconversion  using  two  or  three 
wavelengths.  These  are: 


e.  (>i»  X2)  f'o*  ^g)  4 

( X 1  ,  X  ^ ^  g «  X  2 ) 


M  •  A3'  {  2’  2 

Eneray  and  momentum  conservation  suggest  tnat: 

Tuninc  Mode  > 
Temp . 


9  • 


(X1 ,  X2)  -  (x2,  X] ) 


x;  x. 


h. 


( X 1  ,  X3)  -  (2  x2) 


or  4 
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In  the  diaorams  to  the  rignt  of  e  and  f_  the  dotted  circles  are  the  incut 
Diane  wavelenoths  and  positions  at  the  crossed  circles  represent  the  image 
plane  correlations.  In  e_  the  input  plane  dots  ar,_  crosses  can  be 
interchanged  as  the  alternate  version.  In  f  the  two  inout/outout  plane 
conf iourations  are  shown.  Tne  difference  oetween  temDerature  and  tuninc 
is  that  the  X2  s  are  at  positions  1/3  and  2/3  the  X-j  x^  separation  for 
temperature  tuning  but  are  multiplied  by  constants  (k)  and  (1-k) 
respectively  for  angle  tuning.  Asymmetrical  conversion  (X-j.x^)  -*•(  ,) 
which  wasn't  discussed  above  is  not  the  inverse  of  f.  Similarly, 

( X 2 »  X^)  -*■  ( X -j ,  2X2)  are  more  complex  processes  which  will  not  be  discussed 
here. 
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Section  6 

CONCLUSIONS  AND  RECOMMENDATION 

6. 1  Concl us  ions 

In  this  study  we  have  demonstrated  that  upconversion  correlations  are  being 
controlled  by  the  diffraction  limited  resolution  of  the  lens,  11  ip/mm. 

Since  the  basic  equipment  is  working  well,  the  input  plane  Fourier  transform 
lens  should  be  replaced  using  a  high  quality  50mm  lens  with  an  f/4  aperture 
or  wider.  The  equipment  will  be  ready  then  to  find  the  maximum  resolution 
and  space  bandwidth  product  by  a  variety  of  correlation  methods.  In  order 
of  simplicity  for  integration  into  real  system  applications.  The  correlation 
methods  are: 

1.  Upconversion  -  rotated  input  format 

1.1  (  X-| »  X-j )  -*•  ( ^) 

1.2  ( A-j ,  )  -*•  (  A^) 

2.  Down  Conversion  -  normal  input  format 

2.1  (X-| ,  X2)  -»•  (Aj  ,  A£) 

2.2  (A-j ,  A^)  ■*  (2a2) 

2.3  (Ap  A-j)  -*■  (DC)  Laser  Readout  of  DC  Field 

Joint  transform  correlator  work  has  shown  that  prons  and  liquid  crystal  devices 

have  resolutions  up  to  20  Lp/mm.  With  1  inch  aperture  this  would  correspond 

2 

to  a  SBWP  =  (20  x  25)  .  To  equal  or  exceed  this  value  the  nonlinear  crystal 
would  require  an  input  resolution  of  50  tp/rrm  and  1  cm  aperture.  The  tests 
to  date  suggest  this  is  realistic. 
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The  N-Pattern  test  is  a  very  useful  experimental  tool  which  can  be  developed 
for  optimizing  system  operation. 

6.2  Recommendation 

6.2.1  Finish  Single  Frequency  Upconversion  Experiment  with  Complex  Conjugates 

6. 2. 1.1  Use  NRL  Bit  Shifted  Data 

6. 2. 1.2  Experimental/Theoretical  Study:  2  Frequencies  Up/Down  Conversions 

Both  temperature  and  angle  tuning  (modified  N-test). 

Resolution  vs.  angular  width  of  each  image. 

Resolution  vs.  intersection  angle  2  images. 

6.2.2  Additional  Topics  to  Consider 

6.2.2. 1  Ruby  Laser  Readout  of  "zero  Freq"  Correlation  Field 

6. 2. 2. 2  Two  Crystal  Downcon vers ion  Correlation:  Image  Scaling 

6. 2. 2. 3  For  Completeness:  Reproduce  Russian  3  Crystal  Up/Down  Conversion 
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DEVELOP  EQUATIONS  TO  CHARACTERIZE  FIELD  ANGLE. 


HISTORICAL  SUT'IARY 
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DOWN-CONVERSION  AND  THREE  FREQUENCY  DOWN-CONVERSION. 


FIGURE  1 

APPLICATION  OF  NON-LINEAR  CRYSTALS 
PROGRAM  SUMMARY 
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FIGURE  2 

STANDARD  PATTERNS  OF  CORRELATIONS  AND  CONVOLUTIONS  USING 
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SCHEMATIC  DIAGRAM  OF  CONVOLUTION  AND  CORRELATION  USING  NONLINEAR  OPTICS 


CONVOLUTION  DEFINITIONS 


auto 


FILM-BASED  CORRELATION  OF  PSEUDO-RAIMDOM  CODES 


FILM-BASED  CORRELATION  OF  PSEUDO-RANDOM  CODES 


ti'-fs.'y 


j  'vi*  >■ 
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:  i * 

!  rara 


ZERO-SHIFT  BIT-SHIFT  SHOWING  CORRELATION 

PEAKS  FOR  SECTIONS  MATCHED  TO 
REFERENCE. 


Nd-.YAG  LASER  AND  INPUT  OPTICS  FOR  UP/DOWN  CONVERSION 


RESOLUTION  TESTS  USING  NRL 
35  m  TRANSPARENCIES  OF  BIT  SHIFTED  CODES 


FIGURE  8 

RESOLUTION  TESTS  USI NG  NRL  35  MM  TRANSPARENCIES  OF  BIT  SHIFTED  CODES 

WITH  CD*A  CRYSTAL 


o 


1/100  sec.  e.  1/100  sec. 


CONCLUDE  NO  AUTO-X-CONVOLUTION  SIGNAL  PRESENT 


FIGURE  9 

ROTATION  AND  MASKING  OF  NRL  TRANSPARENCY:  CONVOLUTION  TESTS 


CONVOLUTIONS  WITH  SQUARE  APERTURES 
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In  £  left  most  satellite  has  8  lobe  with  central  blank. 
Others  In  £  may  or  may  not  be  similar. 


CONVOLUTIONS  WITH  SQUARE  APERTURES 


EFFECT  OF  FI  LTERS  ON  TWO  CRYS 


FIGURE  11 

FILTERS  ON  TWO  CRYSTAL  EXPERIMENTS 


Low  Speed  Film 


INTERNAL  REFLECTIONS  AND  S II G 
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Non-parallel  crystal  faces. 

Index  matching  fluid. 

Absorption  filters  in  SHG  Housing. 


FIGURE  12 

INTERNAL  REFLECTIONS  AND  SHG 


SIMULTANEOUS  AUTO-CONVOLUTIONS,  A4-CONVOLUTiONS  AND  NOISE 


CONVOLUTION  NON-LINEARITY 


FIGURE  .  CONVOLUTIONS. 


VARIATION  IN  A  -  A  N  D  AUTO-CONVOLUTION 


FP  PLANE  CRYSTAL:  TEMPERATURE  EFFECTS 


®  #  % 


-CONVOLUTION  VARIATION  WITH  LONGITUDINAL 


Can  be  inverted  to  model  effects  of  crystal  damage  on 
system  degradation. 


N-TEST  OF  CD* A  AT  25°  C 
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NOISE  RELATED  TO  KNOWN  CRYSTAL  DAMAGE 


FIGURE  18 


N-TEST  OF  CD*  A  AT  50°C 


approximation. 


FIGURE  19 
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N-TEST  OF  CD*A  AT  75°C 
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N-TEST  OF  CD* A  AT  100°C 


FIGURE  21 


FIGURE  22 

ANGLE  AND  TRANSLATION  EFFECTS  ON  CD*A  SHG 
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AUTO  -K-  CONVOLUTIONS  WITH  CD*A  CRYSTAL 


FIGURE  23 

AUTO-  X  -CONVOLUTIONS  WITH  CD*A  AT  25°C 
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CONCLUSIONS 


READY  TO  CONSIDER  RUBY  LASER  READOUT  OF  DC  FIELD 


APPENDIX  A-II 
THEORY 


A-II.l  The  Convolution  and  Correlation  Operations 

The  convolution  and  correlation  of  two  input  functions  can  easily  be  formed 
in  an  optical  system.  This  is  done  by  using  the  Fourier  transform  properties 
of  lenses  and  the  square  law  effect  of  an  intensity  sensitive  material  or  a 
nonlinear  device.  The  correlation  function  is  of  particular  interest  since 
it  is  used  in  many  image  identification  and  related  data  processing  operations. 

The  convolution  of  the  function  f(x)  with  g(x)  is  given  by 

oo 

h ( x )  =  f  x  g  =  r  f (u )  g(x  -  u)  du.  (1-1) 

«co 

One  way  of  interpreting  this  is  to  envision  the  function  g  being  reversed  and 
slid  along  f,  as  shown  in  Figure  A-l.  The  convolution  h(x)  is  then  found  to 
be  the  product  of  f  with  g  when  g  is  centered  at  x.  This  product  is  the 
shaded  area  shown  in  the  diagram.  The  new  weighted  distribution  is  a  smoothed 
version  of  the  original. 

One  exercise  that  is  useful  as  a  mental  check  on  the  calculated  convolution  is 
summarized  in  Figure  A-2.  Imagine  that  one  of  the  functions  is  plotted  back¬ 
wards  on  a  movable  piece  of  paper.  For  the  exponential  functions  used  here 
the  product  will  be  zero  when  the  paper  is  to  the  left  of  the  position  shown. 
The  product  slope  then  suddenly  rises  to  some  maximum  value,  where  the  rate 
of  overlap  of  the  two  functions  is  large.  As  the  reversed  function  moves  to 
the  right  of  this  position,  the  convolution  rises  to  a  maximum  and  it  finally 
falls  off  to  zero. 

The  correlation  function  differs  in  that  the  second  function  is  not  reversed 
when  calculating  the  product  area.  The  two  operations  are  the  same,  then,  only 
if  the  input  functions  are  synmetric.  Unlike  the  convolution,  the  correlation 
function  often  results  in  a  peaked  distribution.  This  is  especially  true  when 
one  function  is  similar  to  the  other.  This  peaked  portion  is  due  to  the  large 
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Figure  A-l.  The  convolution  integral  h(x) 
a  shaded  area. 


=  f(x)  *  g(x)  represented  by 


f 

* 

/•« 

•  I 

i! 

L  1 

ryi  i 

i  i 
i  i 

Figure  A-2.  Graphical  construction  for  convolution.  The  movable  piece  of 
paper  has  a  graph  of  one  of  the  functions  plotted  backward. 


area  of  overlap  that  occurs  when  the  functions  are  similar.  The  cross¬ 
correlation  function  is  given  by  the  expression 

f  *  g  ■ 

00 

J ~  f(u)  g*(u  -  x)  du 

= 

00 

J"  f*(u)  g(u  +  x)  du. 

0-2) 

This  function,  being  a  measure  of  similarity  is  a  maximum  when  f  = 
is  the  "autocorrelation"  of  f  with  itself). 

g  (this 

The  formation  of  the  convolution  and  correlation  products  are  also 
Fourier  transforming  operations  by  making  use  of  the  relationships 
Capital  letters  are  used  here  to  denote  the  Fourier  transform  of  a 
function  ( i . e. ,  6  =  F  jg  j ) . 

found  by 

below. 

given 

convolution: 

g  ®  h 

*  F  j G(u)  H(u) |  =  F  Jg(u)  H*(-u)j 

0-3) 

correlation: 

g  *  h  =  F  j G ( u )  H*(u)j 

0-4) 

Because  the  pseudo  random  codes  studied  in  this  contract  are  not 
symmetric,  the  correlation  and  convolution  operations  will  not  be  the  same. 

In  order  to  determine  bit  shift  then,  a  difference  frequency  correlation 
configuration  must  be  used  if  nonlinear  crystals  such  as  CD*A  are  used. 

To  correlate  two  images  with  these  crystals,  a  joint  transform  correlator  can 
be  used.  Here  the  input  functions  are  positioned  a  focal  length  distance  f 
in  front  of  a  Fourier  transforming  lens.  The  Fourier  domain  lies  a  distance 
f  behind  this  lens.  The  spatial  frequencies  that  exist  in  the  input  function 
will  be  displaced  from  each  other  here.  They  will  lie  at  a  distance 

X  =  fx  Xf  (1-5) 
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in  the  Fourier  plane.  Here  f  is  the  frequency  component  in  question,  x  is 
the  optical  processing  wavelength  and  f  is  the  focal  length  of  the  Fourier 
transforming  lens. 

After  the  transforms  are  mixed  by  the  square  law  effect  of  the  crystal,  the 
product  is  inverse  transformed  by  a  second  lens.  According  to  Equation  (4), 
this  yields  the  correlation. 

For  the  case  where  nonlinear  crystals  are  used  to  mix  the  Fourier  domain 
components,  the  output  will  appear  at  a  different  x  than  used  to  illuminate 
the  input.  The  color  of  the  correlation  peak  will  be  the  difference  frequency 
of  the  input. 

The  spatial  extent  required  of  the  nonlinear  crystal  is  governed  by 
Equation  (1-5)  and  the  dimensions  of  the  input  imagery.  For  example,  assume 
that  50  bit/ cm  density  codes  are  to  be  correlated.  In  order  to  obtain  a  well 
resolved  correlation  peak,  frequency  components  several  times  this  value  must 
be  processed.  For  a  large  gain,  let  us  require  that  six  times  the  fundamental 
frequency  be  processed.  Using  Equation  (5),  the  minimum  crystal  width  is 
found : 


x  =  (6  x  50  bits/cm)  (1.06  x  10"4  cm)  (33.02  cm) 

=  1.04  cm  (1-6) 

Since  CD*A  crystals  are  not  grown  much  larger  than  this  1  cm,  a  maximum  bit 
density  of  50  bits/cm  must  be  used.  To  process  data  beyond  this,  some 
modification  of  the  joint  transform  correlation  configuration  must  be  made. 

By  moving  the  input  plane  along  the  optical  axis,  the  size  of  the  Fourier 
transform  can  be  continuously  varied.  Table  A-l  sunmarizes  the  effects  of 
input  plane  position  on  the  phase  and  scale  terms  that  appear  with  the  actual 
transform.  If  the  input  function  is  positioned  a  distance  d  from  the  output 
plane,  then  the  transform  will  be  scaled  by  d,  instead  of  f.  Since  d<f,  the 
transform  will  be  smaller,  and  the  width  required  of  the  CD*A  crystal  will 
not  be  as  extensive.  A  phase  term  will  also  be  introduced,  but  will  be 
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cancelled  when  the  correlation  product  is  formed.  In  this  modified  system, 
therefore,  the  maximum  bit  density  that  can  be  tolerated  at  the  input  will  be 
determined  by  the  crystals  size  and  resolution. 


Table  A-l 


Source 

position 

Object  position 
(relative  to  lensj 

Output 

plane 

Phase 

term 

Scale  propor¬ 
tional  to 

x 

)  in  front 

1  behind 

— 

] 

x 

2, 

J  behind 

*Ji-2.) 

2  f1 

i 

X 

0 

J  behind 

t 

X 

behind  lens  distance  d 
from  output  plane 

/  behind 

x2fld 

j 

A-II.2  Frequency  Conversion 

In  certain  classes  of  crystals  it  is  possible  to  mix  optical  radiation  in 
order  to  generate  sum  and  difference  frequencies.  When  coherent  light  of 
only  one  wavelength  is  available,  the  input  frequency  is  doubled  and  the 
difference  frequency  term  is  zero.  Thus,  harmonic  generation  can  be  used  to 
extend  the  frequency  range  of  available  laser  sources.  Active  or  passive 
parametric  upconversion  can  also  be  used  for  detecting  an  IR  reflected  or 
emitting  scene.  The  IR  radiation  is  imaged  onto  the  nonlinear  upconverting 
crystal.  A  "pump"  laser  then  floods  the  IR  scene  in  the  crystal.  Interaction 
between  the  IR  and  pump  laser  radiation  produces  the  sum  of  the  two  frequencies 
and  a  visible  image  output  is  produced. 

For  our  application  here  the  nonlinear  dielectric  properties  of  these  crystals 
are  used  to  produce  an  output  distribution  that  is  proportional  to  the  square 
of  the  incident  field.  This  phenomena  allows  us  to  obtain  the  convolution  or 
correlation  of  two  input  images,  when  the  crystals  are  used  within  the  Fourier 
domain  of  the  system. 

Let  us  look  at  this  frequency  conversion  operation  in  more  detail.  When  a 
material  is  subjected  to  a  weak  electromagnetic  field,  the  induced  polariza¬ 
tion  is  linear  inside  the  media.  When  higher  optical  field  strengths  are 
present,  such  as  those  available  from  laser  sources,  higher  order  responses 
are  induced.  Of  particular  interest  is  the  second  order  polarization  in  which 
the  sum  or  difference  frequencies  of  two  incident  fields  at  are 

produced.  The  new  polarization  can  be  written  as 


P.(wl  i“2}  -  g  Y\  d...  {_a,l  ^“2*  “1*  -u2}  E  E 

i  o  •  ^  ijk  U-.  i 

J ,K=1 ,2,3  1 


(2-1) 


Here  g  is  a  degeneracy  constant,  d^-k  a  third  order  tensor  of  nonlinear  optical 
coefficients,  and  Ej,  Ek  the  field  amplitudes. 


As  a  result  of  crystal  symmetry,  many  of  the  components  of  d.jk  will  often  be 
zero  or  equal  to  other  components  of  the  crystal. 
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In  order  for  this  new  optical  field  to  propagate,  it  must  travel  at  the  same 
velocity  as  the  incident  fields.  If  this  condition  is  not  met,  the  electro¬ 
magnetic  radiation  created  by  a  dipole  at  one  point  will  not  couple  to  that 
created  elsewhere  and  destructive  interference  will  occur.  In  general,  because 
of  normal  dispersion  this  propagation  velocity  is  different  at  the  sum  and 
difference  frequency  than  it  was  for  the  incident  fields.  The  efficiency  of 
conversion  is  therefore  a  function  of  the  phase  mismatch 


The  wave  vectors  are  equal  to  =  u^  ( 2tt / x n ) ,  where  u  is  a  unit  propagation 
vector. 

One  technique  used  to  satisfy  the  phase-matching  requirement  Ak  -  0  is  to 
take  advantage  of  the  natural  birefringence  of  anisotropic  crystals.  To  see 
how  this  is  done  consider  an  ellipsoid  with  major  axes  in  the  x,  y  and  z 
directions.  This  "index  ellipsoid"  can  then  be  used  to  find  the  two  allowed 
directions  of  polarization  and  the  indices  for  these  directions.  This  is  done 
as  follows.  Through  the  center  of  the  ellipsoid  a  plane  is  drawn  perpendicular 
to  the  direction  of  propagation.  The  ’ntersection  of  this  plane  and  the 
ellipsoid  is  an  ellipse.  The  two  axes  of  this  ellipse  are  parallel  to  the 
two  directions  of  polarization,  and  the  length  of  each  is  equal  to  twice  the 
refractive  index  in  that  direction. 

If  all  three  major  axes  of  the  index  ellipsoid  are  unequal,  then  the  crystal 
is  called  biaxial,  since  two  optic  axes  can  be  defined.  An  optic  axis  is  the 
wave-normal  direction  in  which  the  refractive  index  is  independent  of  the 
direction  of  polarization  (the  direction  perpendicular  to  a  circular  cross 
section  of  the  ellipsoid). 

Uniaxial  crystals  are  more  commonly  encountered.  The  direction  of  polarization 
perpendicular  to  the  optic  axis  is  known  as  the  ordinary  direction.  The  index 
is  independent  of  the  direction  of  propagation  for  such  an  "ordinary"  wave. 

For  the  other  direction  of  polarization,  the  extraordinary  direction,  the  index 
changes  elliptically  between  the  value  of  the  ordinary  index  n  when  the  wave 
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normal  is  parallel  to  the  optic  axis,  and  the  extraordinary  index  ng,  when  the 
wave  normal  is  perpendicular  to  the  optic  axis.  When  the  wave  normal  is  in  a 
direction  e  to  the  optic  axis  the  extraordinary  index  is  given  by 


ne(e)  = 


ne  no 

(n^  sin  2e  +  ng  2cos2e)^2 


(2-3) 


Using  angle  tuning  or  temperature  tuning  of  these  anisotropic  crystals,  the 
phase  matching  requirement  can  be  met. 


Consider,  for  example,  the  case  of  second  harmonic  generation.  Using  Equation 

O 

(2-2)  it  can  be  seen  that  Ak  «  0  if  n  =  n“.  In  normally  dispersive  materials 
the  index  of  the  ordinary  wave  or  the  extraordinary  wave  along  a  given  direc¬ 
tion  increases  with  u.  This  makes  it  impossible  to  satisfy  the  condition 

p 

n  =  n“  when  they  are  both  extraordinary  or  ordinary.  Con:  ider,  however,  a 

o  o 

crystal  with  n„  <  n“.  An  angle  em  then  exists  at  which  n„  (em)  =  n“;  if  the 
fundamental  beam  (at  w)  is  launched  along  em  as  an  ordinary  ray,  the  second- 
harmonic  beam  will  be  generated  along  the  same  direction  as  an  extraordinary 
ray.  This  situation  is  illustrated  in  Figure  A-3.  The  angle  em  is  found  at 
the  intersection  between  the  index  surface  of  the  ordinary  beam  at  u  (circle) 

p 

to  the  index  surface  of  the  extraordinary  ray.  Thi:  gives  n^  (e).  Using 
Equation  (2-3),  for  negative  uniaxial  crystals  (those  for  which  n'^  <  n“)  em 
is  found  to  be 


sin  20 


ip 


(n 


2w,-2 

o  ' 


(n 


2a),-2 

o  ' 


(2-4) 


It  is  seen  that  the  angle  at  which  the  wave  vectors  propagate  determine  the 
phase  mismatch.  These,  and  other  physical  parameters  of  the  system,  can  be 
used  along  with  Maxwell's  equations  to  solve  for  the  intensity  of  the 
upconverted  wave.  If  depletion  of  the  fundamental  wave  is  neglected  and  Ak 
assumed  constant  throughout  the  crystal,  then  the  conversion  efficiency  is 
given  by 


(2-5) 
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*  (optic)  axis 


Figure  A-3.  Normal  (index)  surfaces  for  the  ordinary  and  extraordinary 

rays  in  a  negative  (n£  <  nQ)  uniaxial  crystal.  If  n^“  <  n“, 

the  condition  nf“(e)  =  n“  is  satisfied  at  e  =  e  .  The 
e  o  m 

eccentricities  shown  are  vastly  exaggerated. 


The  phase  mismatch  Ak  due  to  interference  is  seen  to  reduce  efficiency  by  the 
factor  sine  .  The  efficiency  is  also  found  to  be  proportional  to  the 

intensity  PyA  of  the  fundamental  beam. 

So  far  we  have  seen  that  by  proper  orientation  a  crystal  can  be  positioned 
such  that  the  refractive  index  at  a  frequency  u  is  equal  to  the  refractive 
index  at  2u>.  Second  harmonic  generation  can  thus  be  achieved  at  an  angle  e. 

o 

Because  the  angle  at  which  ny*  (e)  =  nQ  varies  with  wavelength,  this  angle 
will  differ  for  different  incident  frequencies.  There  is  one  disadvantage  to 
this  angle  tuning  technique.  Because  the  beams  involved  do  not  propagate 
along  the  principal  axis  of  the  crystal,  the  ordinary  and  extraordinary  beams 
diverge.  This  limits  the  interaction  length,  and  hence  crystal  efficiency. 

An  alternate  means  of  phase  matching  a  crystal  is  by  temperature  tuning. 

This  can  be  done  since  the  extraordinary  index  is  more  temperature  dependent 
than  the  ordinary  index.  This  90°  phase  matching  is  satisfied  by  orientating 
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the  crystal  at  45°  to  both  the  X  and  Y  direction,  and  perpendicular  to  the 
optic  axis.  The  incident  beam  is  polarized  along  the  X-Y  plane,  the  generated 
frequency  along  Z.  These  90°  phase  matched  and  angle  tuning  orientations  are 
shown  in  Figure  A-4  below. 


Figure  A-4.  Crystallographic  orientations  for  angle  and  temperature 
tuning. 

By  looking  at  the  phase  mismatch  versus  e,  the  acceptance  angle  of  a 
particular  configuration  can  be  computed.  In  general.  At  ~  0  only  for  one 
direction  of  the  incoming  beams.  As  rays  come  in  from  angles  other  than 
this,  at  _ecomes  nonzero.  This  results  in  a  decrease  in  efficiency  which 
in  turn  limits  the  field  angle. 

As  an  example,  consider  a  type  I  interaction  in  which  two  ordinary  rays, 
tp  and  tir  intermix  to  form  the  sum  frequency  t$.  The  angular  deviation  of 
the  signal  in  the  plane  of  the  wave  vectors  (parallel  case)  can  then  be 
found: 
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Likewise,  the  acceptance  angle  in  a  direction  perpendicular  to  the  plane  of 
the  wave  vectors  (perpendicular  case)  can  also  be  found  by  wave  vector  analyses, 
vector  analyses.  In  the  small  angle  approximation  it  is  seen  that 


In  general,  it  can  be  seen  that  an  extensive  analysis  is  required  to  find  the 
field  of  view  for  a  particular  system.  The  phase  match  angle  and  temperature 
are  first  found  for  the  central  beams  of  the  reference  and  signal  beams  to  be 
correlated.  The  indices  of  refraction  for  these  beams  are  then  found.  This 
gives  a  numeric  value  to  the  wave  vectors  k^,  kp  and  k$.  An  alternate 
experimental  procedure  is  to  use  the  N-pattern  Test  (Appendix  A-l ,  Section 
4.10)  to  fix  some  of  the  parameters  and  then  complete  the  calculation 
analytical ly. 

The  maximum  allowable  phase  mismatch  is  assumed  to  be  Ak  =  2n/L.  It  is  noted 

that  a  shorter  crystal  is  preferable  to  increase  the  field  of  view.  A  shorter 

crystal,  however,  decrease  the  interaction  length  and  hence  decreases  efficiency. 

Some  optimum  value  must  be  found.  With  the  parameters  k.  ,  k  ,  k  and  Ak,  the 

i  r  p  s 

field  angle  can  be  found. 


A-II.3  Soviet  Papers 

In  the  work  done  by  MDAG  this  year,  the  convolution  of  input  imagery  was 
obtained.  This  was  done  by  using  the  frequency  doubling  properties  of  non¬ 
linear  crystals.  Theoretically,  the  correlation  can  also  be  obtained.  To 
do  this  a  difference  frequency  mode  of  operation  must  be  used.  Since  Soviet 
papers  have  been  published  which  show  the  generation  of  difference  frequency 
harmonics,  these  results  will  be  used  to  compare  to  our  predictions. 

Figure  (A-5)  shows  the  results  published  by  R.  Eremeeva  (Reference  1).  In 
Figure  (A-5),  the  input  data  is  shown  separated  by  a  distance  "a". 


abed 
Figure  A-5.  Soviet  Data 


The  remaining  figures  show  the  results  of  using  different  combinations  of 
input  illumination.  In  Figure  (A-5b)  a  neodymium  laser  at  X1  and  its  third 
harmonic  ( x 3)  are  used;  in  Figure  (A-5c)  X-j  and  x2  are  incident  on  the  input 
images.  The  difference  frequency  correlation  is  recorded  at  the  output.  In 
Figure  (A-5d)  the  convolution  is  found  by  recording  the  second  harmonic  of  x.j 

To  mathematically  describe  the  input  functions,  let  g^  and  g2  be  written  in 
complex  notation  form.  Hence, 

1  /  ~  jw-i  t  juu  t  \ 

Ei"‘?(9l(X'a>e  )  (3 
l/  -JWlt  jWot\  ' 

+  2'(g2^x+a^e  +  92(x+a)e  ) 

The  phase  factors  are  used  to  describe  what  input  wavelengths  illuminate  the 
two  functions  since  u..  =  2ttC/A^  .  A  Fourier  transform  is  next  taken  of  this 
input.  The  scale  and  shift  properties  of  this  operation  are  taken  into 
account.  At  the  Fourier  plane,  the  distribution 
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ef  =  1 


-  jo)-]  t 


-juajV 


G1(u)e’Juae  '  1  +  Gf(u)e"JUtle 

e-: 


,xi 


h2  ua  -J“2t  /  M  J  a2  Ua 

2\-2T  6  g2 \h> 


(3-2) 


is  formed. 


After  passing  through  the  nonlinear  crystal,  this  field  is  squared.  It  is 
then  inverse  transformed  by  a  final  lens.  This  output  is  found  to  be 


2Eq  =  Rii*6(x,w)  +  R2 1 2*^ ( x ' »w) 
+  ^C-|-j*6(x-2a,a)+2a)^ )) 

+  "2'(^'2 1 2 1  (x '  “  ,ou+2co2 )) 


(3-3) 


/  A1  ' 

Ci2i*6lx-a  ^  a ^i^2 
\  2  / 


/  A, 

+  R-J2  i*<5^x+a  +  —  a,a)+ca2-u>i 


+  conjugate  terms. 


The  letters  R  and  C  stand  for  the  correlation  and  convolution  operations, 
respectively.  The  primed  numbers  are  used  to  take  into  account  the  \^/X2 
scale  factor  that  occurs  when  the  input  imagery  is  illuminated  with 
different  wavelengths.  For  example,  the  second  to  last  term  indicates  that 
the  correlation  of  g-|(x)  is  made  with  g2U-|A2x).  This  term  will  appear  at 
a  spatial  location  x  =  -(a  +  X^/X2a),  and  at  a  temporal  frequency  of  (w,-^). 
The  "6"  notation  is  used  to  convey  this  information.  The  conjugate  terms 
mentioned  in  this  equation  are  those  which  are  obtained  when  the  +,  -,  1, 
and  2  signs  are  replaced  with  -,  +,  2  and  1,  respectively. 


1 
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Legend 

1.06  (IR)  =  * 

. 26 ( UV )  =  © 

.53  (green)  =  @ 

.21 (UV)  =  • 

.34  (purple)  =  X 

(a)  Input  X1  =  1.06y,  X3  =  .24u 


v 

/T\ 

1  /"Tv 

X/  _ 

R12 ' 

cn 
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C12‘ 
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cn 
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R12 ' 

(b)  Input  X1  =  1 . 06y ;  X2  =  .53y 
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(c)  Input  X-j  =  1.06y 


©  C22  ©  c22 

VZ7  tz? - - 

cn  ci2  cn 

Figure  A-6.  Predicted  Output  Data  for  Various  Combinations  of  Input 

Illumination.  The  Enclosed  Data  is  Postulated  to  Have  Been 
Photographed  for  Figure  (4). 


This  equation  can  now  be  compared  to  the  experimental  results.  To  do  this, 
the  different  terms  are  plotted  in  Figure  (A-6).  For  example,  x]  and  x3  are 
used  as  inputs.  At  the  output,  four  convolution  terms  appear  in  green,  and 
a  correlation  function  appears  in  violet.  The  term  "oscillating  at  2u>3" 
will  not  appear  since  this  is  beyond  the  response  of  the  available  detectors. 
The  four  functions  in  green  are  those  that  appear  in  the  Soviet  data  of  Figure 
(A-5b). 
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When  X-j  and  x2  are  used  as  the  processing  wavelengths,  then  the  output  data 
appears  as  in  Figure  (A-5b).  Now  four  convolution  functions  appear  in  the  UV, 
one  in  green  and  a  correlation  term  appears  at  1.06y.  Since  the  Soviet  data 
of  Figure  (A-5c)  contain  four  terms,  it  is  assumed  that  their  data  are  a  photo¬ 
graph  of  the  UV  terms  (only  black  and  white  reprints  are  available  for  our 
analysis). 

If  only  one  wavelength  is  present  at  the  input  (X-|  =  x2),  then  Equation  (3-3) 
reduces  to  terms  at  the  three  spatial  locations  x=0,  +  .  This  result  has 
been  verified  experimentally  in  our  lab  and  also  appears  in  Figure  (A-5d). 

It  is  seen  that  only  convolution  terms  are  present.  The  correlation  term  is 
frozen  in  time  since  the  difference  frequency  (o>i  -  u^)  *  0.  Two-frequency 
illumination  is  therefore  preferred  if  correlation  data  is  of  interest. 


A-II.4  Determination  of  Bit  Shift 

It  is  desirable  to  determine  the  number  of  bits  that  one  coded  sequence  is 
shifted  with  respect  to  a  second.  This  is  accomplished  by  using  the  position 
of  the  correlation  peak  at  the  output  of  the  optical  processor. 

The  codes  to  be  processed  have  been  sent  to  us  from  NRL.  According  to  them, 

the  patterns  used  are  generated  using  a  "random  number  generator  that  is 

essentially  a  31  bit  linear  shift  register  with  feedback  taps  to  yield  a 

maximal  length  pseudo  random  sequence.  Such  a  generator  yields  in  excess  of 

g 

2  x  10  bits  before  repeating."  The  codes  are  recorded  in  raster  format  into 
an  M  x  M  array. 

A  total  of  b  bits  are  then  dropped,  and  this  new  bit  shifted  code  is  recorded 
to  the  right  of  the  reference  code. 

A-II.4.1  Correlation  of  Zero-Bit  Shifted  Imagery 

The  correlation  of  two  bit  images  can  be  predicted  by  decomposing  each  image 
into  the  sum  of  its  individual  bits.  For  example,  let  g  be  the  M  x  M  reference 
array  composed  of  the  bits  g...  Then 

*  J 

M  M 

'  ’EE  <4-» 

i=l  j=l 

Similarily,  let  the  second  image  be  defined  as 

M  M 

h  ■  Z  £  v  (4-2> 

k=l  *«1 

The  correlation  of  the  two  images  is  then  given  by  the  transformed  product 
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(4-3) 


■  F''  {?  J  G,j  "«*  ‘j2,(u4x  *  M 

+  cross  terms 

In  the  latter  expression  ax  and  Ay  are  the  separation,  at  the  input  plane, 
between  the  two  bits  of  g  and  h  of  indices  (i,  j).  Since  the  correlation  of 
these  bits  overlap  at  the  output  plane,  the  other  terms  are  much  lower  in 
intensity  by  comparison.  These  "noise"  terms  are  therefore  neglected. 

After  carrying  out  the  inverse  Fourier  transform  of  Equation  (4-3)  the  location 
of  the  correlation  peak  is  found  to  be 

M  M 

9  *  h  =  £  £  9^  *  h- .  *  6(x  -  ax,  y  -  Ay)  (4-4) 

1=1  j=l  J  IJ 

For  two  M  x  M  arrays  spaced  a  distance  "a"  from  each,  define  the  coordinate 
system  shown  in  Figure  A-6. 

(M  +  a,  1) 


y 


Figure  A-6.  Coordinate  system  for  zero  bit  shifted  array. 

For  this  case  Ax  =  M  +  a.  Ay  =  0.  The  correlation  peak  is  therefore  located 
at  a  relative  output  position  of  (m  +  a,  0), 
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All. 4. 2  Correlation  of  Bit  Shifted  Images 


Consider  next  the  situation  where  a  total  of  b  bits  have  been  dropped  from 
the  right  hand  image  before  It  was  recorded.  In  Figure  (A-7)  the  new 
reference  and  signal  Images  are  depicted.  The  hatched  area  represents  the 
bits  that  have  been  dropped  to  form  the  signal  image.  Areas  numbered  1  and  2 
are  matched. 


0,  1) 

(R(b/M)  +  1,  L (b/M)  +  1) 


(M  +  a,  1) 

(2M  +  a  -  R(b/M)  +1,1) 


(1 ,L(b/m)  +  2) 


HR 

TA 

l 

© 

© 

© 

© 

REFERENCE 

SIGNAL 

(2M  +  a,  M) 


Figure  A-7.  Coordinates  of  1st  Bit  of  Matched  Imagery  Sections. 


Using  Equation  (4-3)  the  location  of  the  correlation  peak  is  seen  to  be  a 
function  of  the  separation  between  matched  data  sections.  To  determine  the 
values  of  ax  and  Ay  for  the  respective  image  sections,  define  the  coordinate 
of  the  first  bit  (of  reference)  to  be  used  in  signal  image  as 

( R ( b/m)  +  1,  L(b/m)  +1)  (4- 

Here  the  symbol  "L"  is  used  to  denote  the  "floor"  of  the  quotient  b  *  m,  or 
the  next  smallest  integer  "R"  is  the  remainder  for  this  product.  Coordinates 
are  then  assigned  to  the  Images. 
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The  correlation  of  the  bit  shifted  image  can  now  be  predicted: 


9  *  h  “  (9©*  h0>  +  (9©  *  h©> 

*  9  Q  *  hQ  *  «(M  +  a  -  R(b/M)-1,  -L(b/M) 

+90*  5 (2M  +  a  -  R(b/M) ,  -L ( b/M )  -1) 

(4-6j 

The  magnitude  of  these  two  terms  will  be  proportional  to  the  areas  of  images 
©  and  (2)  ,  respectively. 

As  an  example,  consider  the  case  of  a  5  x  5  bit  array  where  the  first  7  bits 
have  been  dropped.  The  coordinates  of  the  matched  data  sections  are  found  by 
reference  to  Figure  A-7.  The  initial  coordinates  of  the  first  reference  bit 
found  in  the  signal  is  (R(b/M)  +1,  L ( b/M)  +1)  *  (R(7/5)  +  1,  L(7/5)  +  1)  *  (3,2). 
The  first  coordinate  of  signal  section  No.  2  is  (2M  +  a  -  R(b/M)  +1,1)* 

(2.5  +  a  -  R(7/5)  +  1,  1)  *  (9  +  9,  1).  This  is  shown  in  Figure  A-8. 


Figure  A-8.  Coordinates  for  7  Bits  Shifted  on  5  x  5  Array. 
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From  these  values  the  spatial  shifts  ax  and  Ay  are  found  for  the  respective 
image  sections.  This  shifts  the  correlation  peak  into  the  two  positions  given 
by  Equation  (4-5).  In  this  case  g^  *  h^  is  found  at  R^(2  +  a,  -1),  and 
92  *  h 2  at  ^2^  +  Tbe  locations  of  these  two  correlation  points  are 

also  plotted  in  Figure  A-8. 

For  the  film-based  correlations  shown  in  Figure  4,  Appendix  A-I,  a  wrap¬ 
around  format  was  used  to  form  the  bit-shifted  signal.  For  this  reason 
four  images  matched  up  with  reference  and  four  corresponding  correlation 
peaks  were  formed. 
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